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152. 


A MEMOIR ON THE THEORY OF MATRICES. 


[From the Philosophical Transactions of the Royal Society of London, vol. cxivit. for 
the year, 1858, pp. 17—37. Received December 10, 1857,—Read January 14, 1858.] 


THE term matrix might be used in a more general sense, but in the present 
memoir I consider only square and rectangular matrices, and the term matrix used 
without qualification is to be understood as meaning a square matrix; in this restricted 
sense, a set of quantities arranged in the form of a square, e.g. 


(64 04 © 2 
Gi Oded 
a, 05 o 


is said’ to be a matrix. The notion of such a matrix arises naturally from an 
abbreviated notation for a set of linear equations, viz. the equations 
X=au +by +cz, 
Y=aa +by +02, 
Z = aw + by + ¢”z, 
may be more simply represented by 
(X, Y, Z)=(a, b, ¢ Qa, y, 2), 
a’, b’, € 
a” b” o 
and the consideration of such a system of equations leads to most of the fundamental 
notions in the theory of matrices. It will be seen that matrices (attending only to 


those of the same order) comport themselves as single quantities; they may be added, 


60—2 
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multiplied or compounded together, &c.: the law of the addition of matrices is pre- 
cisely similar to that for the addition of ordinary algebraical quantities; as regards 
their multiplication (or composition), there is the peculiarity that matrices are not in 
general convertible; it is nevertheless possible to form the powers (positive or negative, 
integral or fractional) of a matrix, and thence to arrive at the notion of a rational 
and integral function, or generally of any algebraical function, of a matrix. I obtain 
the remarkable theorem that any matrix whatever satisfies an algebraical equation of 
its own order, the coefficient of the highest power being unity, and those of the 
other powers functions of the terms of the matrix, the last coefficient being in fact 
the determinant; the rule for the formation of this equation may be stated in the 
following condensed form, which will be intelligible after a perusal of the memoir, 
viz. the determinant, formed out of the matrix diminished by the matrix considered 
as a single quantity involving the matrix unity, will be equal to zero. The theorem 
shows that every rational and integral function (or indeed every rational function) of 
a matrix may be considered as a rational and integral function, the degree of which 
is at most equal to that of the matrix, less unity; it even shows that in a sense, 
the same is true with respect to any algebraical function whatever of a matrix. One 
of the applications of the theorem is the finding of the general expression of the 
matrices which are convertible with a given matrix. The theory of rectangular 
matrices appears much less important than that of square matrices, and I have not 
entered into it further than by showing how some of the notions applicable to these 
may be extended to rectangular matrices. 


1. For conciseness, the matrices written down at full length will in general be 
of the order 3, but it is to be understood that the definitions, reasonings, and con- 
clusions apply to matrices of any degree whatever. And when two or more matrices 
are spoken of in connexion with each other, it is always implied (unless the contrary 
is expressed) that the matrices are of the same order. 


2. The notation 
C.o, Bc Ba, ¥,' 2) 


represents the set of linear functions 


((a, b, cha, y, 2), (a, b, cha, y, 2), (a, b”, cYa, y, 2)), 
so that calling these (X, Y, Z), we have 
(X, Y, Z)=(a, b, c Ga, y, 2) 
Big ty, 18 


| a. b”, c” | 


and, as remarked above, this formula leads to most of the fundamental notions in the 
theory. 


www.rcin.org.pl 


152] A MEMOIR ON THE THEORY OF MATRICES. 477 


3. The quantities (X, Y, Z) will be identically zero, if all the terms of the matrix 
are zero, and we may say that 


(0, F00 
OERO 
028 98 


is the matrix zero. 


Again, (X, Y, Z) will be identically equal to (a, y, z), if the matrix is 


Ct, By, 0) 
oe: See 
ee ees | 


and this is said to be the matrix unity. We may of course, when for distinctness it 
is required, say, the matrix zero, or (as the case may be) the matrix unity of such an 
order. The matrix zero may for the most part be represented simply by 0, and the 
matrix unity by 1. 
4, The equations 
(X, ¥, 2=(a, 6, ¢ Ya, y, 2), X,Y, Z=(a, B, y Gey, 2) 


a”, ae Cc” ar. B y” 
give 
(X +X, Y+Y,Z+Z)=(a +a, b+ß8, c +y Qa. y, z) 
a’ +a’, b’ +f’, c +y 
a” +a", b” +B”, o” +y” 
and this leads to 
e A +8),00 #78 aE T 0, 00 E E: E y ) 
a 4. a’ : b a jeg c 4 y a’ y b’ ‘ c a’ : Bi $ y 
a” + a’, br +P", c” +y” a”, b”, g a, BS y” 
as a rule for the addition of matrices; that for their subtraction is of course similar 
to it. , 


5. A matrix is not altered by the addition or subtraction of the matrix zero, 
that is, we have M + 0 = M. 


The equation L= M, which expresses that the matrices L, M are equal, may also 
be written in the form ZL—M=0, i.e. the difference of two equal matrices is the 
matrix zero. 


6. The equation Z=— M, written in the form L+ M=0, expresses that the sum 
of the matrices Z, M is equal to the matrix zero, the matrices so related are said to be 
opposite to each other; in other words, a matrix the terms of which are equal but oppo- 
site in sign to the terms of a given matrix, is said to be opposite to the given matrix. 
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7. It is clear that we have L+ M= M + L, that is, the operation of addition is 
commutative, and moreover that (2+M)+N=L4+(M+N)=L+M+VN, that is, the 
operation of addition is also associative. 


8. The equation 
(X; Y, A=(a, bp c Yama, my, mz) 


uy b” mA 
> 
written under the forms 


(X, Y, 2=m(a, b, ¢ Yay, z)=( ma, mb, mo Ja, y, 2) 


ioe Aes ma’, mb’, me’ 
gives 


ma’, mb’, me 


a”, b v ar K mb”, mec 


as the rule for the multiplication of a matrix by a single quantity. The multiplier m 
may be written either before or after the matrix, and the operation is therefore com- 
mutative. We have it is clear m (L + M)=mL+mM, or the operation is distributive. 


9. The matrices L and mL may be said to be similar to each other; in 
particular, if m=1, they are equal, and if m = — 1; they are opposite. 


10. We have, in particular, 
m(1, 0, O)=(m, 0, 0), 
EUT 3) 4 0, m, 0 
0, 5D, wk 0, 0, m 
or replacing the matrix on the left-hand side by unity, we may write 
m=(m, 0, 0); 
0, m, 0 
0, 0, m 
the matrix on the right-hand side is said to be the single quantity m considered as 


involving the matria unity. 


11. The equations 


(X, Y, D=(a, b, ¢ Wy, 2; @ y za=(a, B, y BE 9, 
Pa ee a, Byoce 
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give 
(X,Y, 2)=(A, B, © En =a, b, c a, B, y Emt) 
Bee: ae Cs OAS hanai Bh 
Ja oe? Krik i ey Oo" -UR OE N 
and thence, substituting for the matrix 
A ya B bichon) 
6: Te o Seem Hs 


A B Ee 

its value, we obtain 
((a,b,c a d,a”), (a,b,c U8, 6,8"), (a,b,c Sy Y y) )=(a, b,c ha, By) 
(a’, b’ ‘ c Ya, a’, a”), (a’ d b’ ; č UB, B’, 8B”), (ars b’ } c ýy, Y, y”) a”, y i g a’ l B í y 
Ck, b”, o” a, a’, a”), la ike c” UR, B’, BY (a”, b”, c” ýy, Y, y”) a bs oe a”, a y” 
as the rule for the multiplication or composition of two matrices. It is to be 
observed, that the operation is not a commutative one; the component matrices may 
be distinguished as the first or further component matrix, and the second or nearer 
component matrix, and the rule of composition is as follows, viz. any line of the com- 
pound matrix is obtained by combining the corresponding line of the first or further 
component matrix successively with the several columns of the second or nearer com- 
pound matrix. 

[We may conveniently write 


(2, a, a"), (B, B, BY), (wy ¥) 


(a ’ b , c ) » 9 3 
(a’ ? b , A ) » » » 
(a”, as c”) ” ” ” 


to denote the left-hand side of the last preceding equation. ] 


12. A matrix compounded, either as first or second component matrix, with the 
matrix zero, gives the matrix zero. The case where any of the terms of the given 
matrix are infinite is of course excluded. 


13. A matrix is not altered by its composition, either as first or second component 
matrix, with the matrix unity. It is compounded either as first or second component 
matrix, with the single quantity m considered as involving the matrix unity, by 
multiplication of all its terms by the quantity m: this is in fact the before-mentioned 
rule for the multiplication of a matrix by a single quantity, which rule is thus seen 
to be a particular case of that for the multiplication of two matrices. 


14. We may in like manner multiply or compound together three or more 
matrices: the order of arrangement of the factors is of course material, and we may 
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distinguish them as the first or furthest, second, third, &c., and last or nearest 
component matrices: any two consecutive factors may be compounded together and 
replaced: by a single matrix, and so on until all the matrices are comipounded together, 
the result being independent of the particular mode in which the composition is 
effected; that is, we have L.MN=IM.N=LMN, LM.NP=L.MN.P, &c., or the 
operation of multiplication, although, as already remarked, not commutative, is associative. 


15. We thus arrive at the notion of a positive and integer power L? of a 
matrix L, and it is to be observed that the different powers of the same matrix are con- 
vertible. It is clear also that p and q being positive integers, we have L?. [9= [?+4, 
which is the theorem of indices for positive integer powers of a matrix. 


16. The last-mentioned equation, Z? . L3 = +9, assumed to be true for all values 
whatever of the indices p and g, leads to the notion of the powers of a matrix for any 
form whatever of the index. In particular, Z’. L° = L? or L’ =1, that is, the Oth power 
of a matrix is the matrix unity. And then putting p=1, g=—1, or p=—1, g=1, we 
have L.L~=L>.L=1; that is, J~, or as it may be termed the inverse or reciprocal 
matrix, is a matrix which, compounded either as first or second component matrix 
with the original matrix, gives the matrix unity. 


17. We may arrive at the notion of the inverse or reciprocal matrix, directly 
from the equation 


(X, ¥, A=(a, b, c Re, y, 2), 


a’ 2 b ; c’ 


in fact this equation gives 
(æ, y, z)=( A, AY, A” UX, Y, Z=((a, b, ¢ RY, Y, J), 
| 
Bey ae 
Gy sf 


Kio © 


a’ y b”, c” 


and we have, for the determination of the coefficients of the inverse or reciprocal 
matrix, the equations 


(A, A, A” Ya, by o )=(1, 0, 0), 


B, B, Bila. ve a 

C, O g g” y b”, e” f 0, $ 

E T EART, ea Bogs MERA Epa 
K b, ¢ || B, B, B’ Oe ERT: 
Large cel ee ees eee 
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which are equivalent to each other, and either of them is by itself sufficient for the 
complete determination of the inverse or reciprocal matrix. It is well known that if 


V denote the determinant, that is, if 


Vela, ib, o 


a b”, e” 


then the terms of the inverse or reciprocal matrix are given by the equations 


Hik, 1 CORIO PRA OS O11, ee 
Vv 0 Oa eo x A ART 
0, b”, o” a”, 0, o” 


E E O E n E a pl 
D A i V Tae Bev. N 
oe, 8. OT OV, OV, OV 


where of course the differentiations must in every case be performed as if the terms 
a, b, &c. were all of them independent arbitrary quantities. 


18, The formula shows, what is indeed clear @ priori, that the notion of the 
inverse or reciprocal matrix fails altogether when the determinant vanishes: the matrix 
is in this case said to be indeterminate, and it must be understood that in the 
absence of express mention, the particular case in question is frequently excluded from 
consideration. It may be added that the matrix zero is indeterminate; and that the 
product of two matrices may be zero, without either of the factors being zero, if only 
the matrices are one or both of them indeterminate. 


19. The notion of the inverse or reciprocal matrix once established, the other 
negative integer powers of the original matrix are positive integer powers of the 
inverse or reciprocal matrix, and the theory of such negative integer powers may be 
taken to be known. The theory of the fractional powers of a matrix will be further 
discussed in the sequel. 


20. The positive integer power L™ of the matrix L may of course be multiplied 
by any matrix of the same degree: such multiplier, however, is not in general. con- 
vertible with Z; and to preserve as far as possible the analogy with ordinary 
algebraical functions, we may restrict the attention to the case where the multiplier 
is a single quantity, and such convertibility consequently exists. We have in this 
manner a matrix cL™, and by the addition of any number of such terms we obtain 
a rational and integral function of the matrix L. 

C. II 61 


www.rcin.org.pl 


482 A MEMOIR ON THE THEORY OF MATRICES. [152 


21. The general theorem before referred to will be best understood by a com- 
plete development of a particular case. Imagine a matrix 


M=(a, 6), 
c, a 
and form the determinant 
a— M, b : 
c ,a—-M 


the developed expression of this determinant is 
Me — (a + d) M + (ad — be) M’; 
the values of M?, M, M° are 
(.at+ G n bda+d).), rf ar arD, a S); 
|c@+dp, +o | |e, al [o 1| 


and substituting these values the determinant becomes equal to the matrix zero, viz. 
we have 


a— M, b =(a@+be , b(at+d) )-(a+d)(a, b)+(ad—be) (1, 0) 
eee oh gate Bebe | |o al lo, 1| 
=( (a+ be)—(a+d)a+(ad—be), b(a+d)—(a+d)b y=(_0,~ 0); 
E MP he-(a4d)d+ad—bo| O ef 
that is 
a— M, b = 0), 
C , a—-M 


where the matrix of the determinant is 
( a, b )- M ( l, 0 ), 
Pacala he 


that is, it is the original matrix, diminished by the same matrix considered as a single 
quantity involving the matrix unity. And this is the general theorem, viz. the deter- 
minant, having for its matrix a given matrix less the same matrix considered as a 
single quantity involving the matrix unity, is equal to zero. 


22. The following symbolical representation of the theorem is, I think, worth 
noticing: let the matrix M, considered as a single quantity, be represented by M, then 
writing 1 to denote the matrix unity, M.1 will represent the matrix M, considered 
as a single quantity involving the matrix-unity. Upon the like principles of notation, 
I.M will represent, or may be considered as representing, simply the matrix M, and 


the theorem is A nb 
Det. (1. M-—M.1)=0. 
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23. I have verified the theorem, in the next simplest case of a matrix of the 
order 8, viz. if M be such a matrix, suppose 


M=(a, b, c), 


dye 6,~ fF 
gio h, | 


then the derived determinant vanishes, or we have 


a—M, b jek = 0, 
d , e-M, f 
g Rn , t-M 


or expanding 
M -—(a+e+i) M+ (ei +a +ae— fh —cg — bd) M — (aei + bfg + cdh — afh — bdi — ceg) =0; 


but I have not thought it necessary to undertake the labour of a formal proof of 
the theorem in the general case of a matrix of any degree. 


24. If we attend only to the general form of the result, we see that any matrix 
whatever satisfies an algebraical equation of its own order, which is in many cases the 
material part of the theorem. 


25. It follows at once that every rational and integral function, or indeed every 
rational function of a matrix, can be expressed as a rational and integral function of 
an order at most equal to that of the matrix, less unity. But it is important to 
consider how far or in what sense the like theorem is true with respect to irrational 
finctions of a matrix. If we had only the equation satisfied by the matrix itself, 
such extension could not be made; but we have besides the equation of the same 
order satisfied by the irrational function of the matrix, and by means of these two 
equations, and the equation by which the irrational function of the matrix is deter- 
mined, we may express the irrational function as a rational and integral function of 
the matrix, of an order equal at most to that of the matrix, less unity; such expression 
will however involve the coefficients of the equation satisfied by the irrational function, 
which are functions (in number equal to the order of the matrix) of the terms, 
assumed to be unknown, of the irrational function itself. The transformation is never- 
theless an important one, as reducing the number of unknown quantities from n? (if n 
be the order of the matrix) down to n. To complete the solution, it is necessary to 
compare the value obtained as above, with the assumed value of the irrational function, 
which will lead to equations for the determination of the n unknown quantities. 


26. As an illustration, consider the given matrix 


M=(a, b), 


le, a| 
61—2 
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and let it be required to find the matrix L=VM. In this case M satisfies the equation 


M—(a+d)M+ad—be=0; 
and in like manner if 
L=(a B) 
y, È 
then Z satisfies the equation 
I? — (a + è) L+ að- By=0; 


and from these two equations, and the rationalized equation Z? =M, it should be possible 
to express Z in the form of a linear function of M: in fact, putting in the last 
equation for L? its value (= M), we find at once 


1 
L= rai [M + (aò — By)], 
which is the required expression, involving as it should do the coefficients 2+ 8, ad— Sy 


of the equation in ZL. , There is no difficulty in completing the solution; write for 
shortness a+ =X, aô— y= Y, then we have 


L=(a B)=(at+Y b ), 
uf 3 | r X 

KA d+Y 
P AESA y 


and consequently forming the values of a +ò and aò — By, 


y=2%+4+2Y 
X , 
ya @tV a+ Y)—be 
aeee ad” Rote! vd 


and putting also a+ d= P, ad—bc=Q, we find without difficulty 
X=J/P+29Q, 
Y=VQ, 


and the values of a, 8, y, 8 are consequently known. The sign of VQ is the same in 
both formule, and there are consequently in all four solutions, that is, the radical VM 
has four values, 


27. To illustrate this further, suppose that instead of M we have the matrix 


M?=(a, b Y=( @+be, b(a+d) ), 
P d | | o(a+d), d? + be 
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so that L? = M?, we find 
P=(a+dy—2(ad — be), 
Q = (ad — bey, 
and thence V Q=+ (ad —bc). Taking the positive sign, we have 
Y= ad-—be, 
X=+(a+d), 
and these values give simply 
L=+( a, b )=4M. © 
lo d | 


But taking the negative sign, 
Y =—ad + bo, 


X =+ V (a — d? + 4be, 
and retaining X to denote this radical, we find 


pat a? — ad + 2bc b(a+d) ), 


x x 
c(a +d) d? — ad + 2bc 
a j X 


which may also be written 
-2+9 (a, b) 2(ad—be) (1, 0), 
X |o a | sod oe A MA 


or, what is the same thing, 
Latte y_ 2a — be), 
x r Sa 


and it is easy to verify à posteriori that this value in fact gives Z27=M% It may 
be remarked that if 


M=( 1, 0 }=1, 
0, 1| 
the last-mentioned formula fails, for we have X=0; it will be seen presently that 


the equation Z?=1 admits of other solutions besides L=+1. The example shows how 
the values of the fractional powers of a matrix are to be investigated. 


28. There is an apparent difficulty connected with the equation satisfied by a 
matrix, which it is proper to explain. Suppose, as before, 


M =( a, b ), 
lo d | 
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so that M satisfies the equation 
| a—M, b =0, 


le , a-M 
M —(a+d)M+ad—be=0, 


or 


and let X, X, be the single quantities, roots of the equation 


a— x, b =0 
c , a-X 


or 


X?—(a +d) X +ad—be =0. 
The equation satisfied by the matrix may be written 
(M = X, (M= xX,)= 0, 


in which X, X, are to be considered as respectively involving the matrix unity, and it 
would at first sight seem that we ought to have one of the simple factors equal to 
zero; this is obviously not the case, for such equation would signify that the perfectly 
indeterminate matrix M was equal to a single quantity, considered as involving the 
matrix unity. The explanation is that each of the simple factors is an indeterminate 
matrix, in fact M — X, stands for the matrix 


ho 0, ah 
0) Wh eae a 


and the determinant of this matrix is equal to zero. The product of the two factors 
is thus equal to zero without either of the factors being equal to zero. 


29. A matrix satisfies, we have seen, an equation of its own order, involving the 
coefficients of the matrix; assume that the matrix is to be determined to satisfy some 
other equation, the coefficients of which are given single quantities. It would at first 
sight appear that we might eliminate the matrix between the two equations, and thus 
obtain an equation which would be the only condition to be satisfied by the terms 
of the matrix; this is obviously wrong, for more conditions must be requisite, and we 
see that if we were then to proceed to complete the solution by finding the value of 
the matrix common to the two equations, we should find the matrix equal in every case 
to a single quantity considered as involving the matrix unity, which it is clear ought 
not to be the case. The explanation is similar to that of the difficulty before adverted 
to; the equations may contain one, and only one, common factor, and may be both of 
them satisfied, and yet the common factor may not vanish. The necessary condition 
seems to be, that the one equation should be a factor of the other; in the case where 
the assumed equation is of an order equal or superior to the matrix, then if this 
equation contain as a factor the equation which is always satisfied by the matrix, the 
assumed equation will be satisfied identically, and the condition is sufficient as well 
as necessary: in the other case, where the assumed equation is of an order inferior 
to that of the matrix, the condition is necessary, but it is not sufficient. 
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30. The equation satisfied by the matrix may be of the form M"™=1; the 
matrix is in this case said to be periodic of the nth order. The preceding conside- 
rations apply to the theory of periodic matrices; thus, for instance, suppose it is 
required to find a matrix of the order 2, which is periodic of the second order. Writing 


M=(a, 6), 
ead 
we have 
M? -— (a+d)M + ad-— be =0, 
and the assumed equation is 
M?—1=0. 

These equations will be identical if 

a+d=0, ad—be=-—1, 


that is, these conditions being satisfied, the equation M?—1=0 required to be satisfied, 
will be identical with the equation which is always satisfied, and will therefore itself 
be satisfied. And in like manner the matrix M of the order 2 will satisfy the 
condition M*—1=0, or will be periodic of the third order, if only M*—1 contains as 
a factor 

M? —(a+d)M+ ad — be, 


and so on. 


31. But suppose it is required to find a matrix of the order 3, 


M=(-a,. b, ¢) 
AS 
yh, NE 
which shall be periodic of the second order. Writing for shortness 
a-M, 6b ,c  |=-(f-—AlP+BM-O, 
d , e-M, f 
g w ROY M 


the matrix here satisfies 
M:—AM?+ BM —C=0, 


and, as before, the assumed equation is M?—1=0. Here, if we have 1+B=0, 4+C=0, 
the left-hand side will contain the factor (M?— 1), and the equation will take the form 
(M?-—1)(M+C)=0, and we should have then M?—1=0, provided M+C were not an 
indeterminate matrix. But M+ C denotes the matrix 


(a+, b > c ) 
e oy een y 
g 3, Vb OS ae 
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the determinant of which is C*+AC?+BC+C, which is equal to zero in virtue of 
the equations 1+ B=0, A+C=0, and we cannot, therefore, from the equation 
(M*?—1)(M+C)=0, deduce the equation M?-1=0. This is as it should be, for the 
two conditions are not sufficient, in fact the equation 
M*=(a?+bd+cg, ab+be+ch, ac+bf+ci )=1 
da+ed+fg, db+e+fh, de+ef+fr 
gathd+ig, gbt+the+th, got+thf+? 


gives nine equations, which are however satisfied by the following values, involving in 
reality four arbitrary coefficients; viz. the value of the matrix is 


Sear (Btu: (B+) ee 
atte.” at+B+y ° at+B+y 
— (y + a) wy B —(y +a) Aw 

St Behe] 5 Ce Ry | a+B+y 
—(@+B)my* = —(@+ 8) vA> y 
at+B+y ’ Rb BRM ST nO Bae ig 


so that there are in all five relations (and not only two) between the coefficients of 
the matrix. 


32. Instead of the equation M"—1=0, which belongs to a periodic matrix, it is 
in many cases more convenient, and it is much the same thing to consider an 
equation M"—k=0, where k is.a single quantity. The matrix may in this case be 
said to be periodic to a factor près. s 


33. Two matrices L, M are convertible when LM = ML. If the matrix M is given, 
this equality affords a set of linear equations between the coefficients of L equal in 
number to these coefficients, but these equations cannot be all independent, for it is 
clear that if Z be any rational and integral function of M (the coefficients being single 
quantities), then Z will be convertible with M; or what is apparently (but only appa- 
rently) more general, if L be any algebraical function whatever of M (the coefficients 
being always single quantities), then Z will be convertible with M. But whatever the 
form of the function is, it may be reduced to a rational and integral function of an 
order equal to that of M, less unity, and we have thus the general expression for 
the matrices convertible with a given matrix, viz. any such matrix is a rational and 
integral function (the coefficients being single quantities) of the given matrix, the 
order being that of the given matrix, less unity. In particular, the general form of 
the matrix LZ convertible with a given matrix M of the order 2, is L=aM+8, or 
what is the same thing, the matrices 


(a, 6), (a, 0) 
| 6, ah! le, d 


will be convertible if -d :0':c=a-—d:b:c. 
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34. Two matrices L, M are skew convertible when LM =— ML; this is a relation 
much less important than ordinary convertibility, for it is to be noticed that we cannot 
in general find a matrix L skew convertible with a given matrix M. In fact, con- 
sidering M as given, the equality affords a set of linear equations between the coeffi- 
cients of L equal in number to these coefficients; and in this case the equations are 
independent, and we may eliminate all the coefficients of L, and we thus arrive at a 
relation which must be satisfied by the coefficients of the given matrix M. Thus, 
suppose the matrices 


(a5 6), Cah 68) 
PIPN ai a’ 


are skew convertible, we have 


Sa Bie he aa He ek oes 
Le AY cool tee a 


(a, 0 )( a, b )=(aad'+bc, wb+dd ), 
| C. 7a | | c, d | | cat+d'c, cb +d'd | 
and the conditions of skew convertibility are 
2aa’ + be’ + b’c =0, 
b'(a +d) +b(a' +d’)=0, 
c (a +d) +c(a’' +d’) =0, 
2dd’ + be’ + b’c = 0, 


Eliminating a’, b, c’, d’, the relation between a, b, c, d is 


2a, ee O is Air 
b, a+d, ; b 
Cc, a+d, c 

go. I Bd 


which is 
(a + dY} (ad — be) = 0. 


Excluding from consideration the case ad—be=0, which would imply that the matrix 

was indeterminate, we have a+d=0. The resulting system of conditions then is 
a+d=0, a +d'=0, aa +be'+b'c+dd' =0, 

the first two of which imply that the matrices are respectively periodic of the second 


order to a factor près. 


35. It may be noticed that if the compound matrices LM and MLZ are similar, 
they are either equal or else opposite; that is, the matrices Z, M are either convertible 
or skew convertible. 


o AD i a 62 
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36. Two matrices such as 
CaO. eerie Ty 
ERE de Maa st ae 


are said to be formed one from the other by transposition, and this may be denoted 
by the symbol tr.; thus we may write 


(-a,. c ):=tr..(.a, 0_). 
ls a! |e al 
The effect of two successive transpositions is of course to reproduce the original matrix. 
37. It is easy to see that if M be any matrix, then 
(tr. MX = tr. (MM), 
(tr. M)> = tr. (M—). 


and in particular, 


38. If L, M be any two matrices, 
tr. (LM) =tr. M. tr. L, 
and similarly for three or more matrices, L, M, N, &c., 


tr. (ZMN) =tr. N. tr. M. tr. L, &e. 


40. A matrix such as 


C4; JA "gd 
l a 
gi fee 


which is not altered by transposition, is said to be symmetrical. \ 
41. A matrix such as 
PFD vy, a, 
— yp, 0, ` 
H —X, 0 


which by transposition is changed into its opposite, is said to be skew symmetrical. 


42. It is easy to see that any matrix whatever may be expressed as the sum of 
a symmetrical matrix, and a skew symmetrical matrix; thus the form 


( a, h+y», g—-p ) 
h-v, b, fN 
Gti ba o 8 


which may obviously represent any matrix whatever of the order 3, is the sum of the 
two matrices last before mentioned. 
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43. The following formulæ, although little more than examples of the composition 
of transposed matrices, may be noticed, viz. 
a, bYa, c)=( AFO, ac+hd ) 
d, b | | ac+bd, E+ | 
which shows that a matrix compounded with the transposed matrix gives rise to a 
symmetrical matrix. It does not however follow, nor is it the fact, that the matrix and 
transposed matrix are convertible. And also 
(a, c{a ba c )=(a@+bed+a(h+e), oc +abd+c(a?+d ) 
Le alle, N E s-e+ad+e@t+d), @+abo+dh+e) 


which is a remarkably symmetrical form. It is needless to proceed further, since it 


is clear that 
(a, c¥a bYa cha b)=((a, cha, b Yy: 
FNP ae ae ee se ee et 
44. In all that precedes, the matrix of the order 2 has frequently been con- 
sidered, but chiefly by way of illustration of the general theory; but it is worth while to 


develope more particularly the theory of such matrix. I call to mind the fundamental 
properties which have been obtained, viz. it was shown that the matrix 


M (sa, 0.2), 
le, a | 


satisfies the equation 
M*?—(a+d)M+ad—be=0, 


and that the two matrices 


will be convertible if 
a—d:b':c=a-d:b:¢, 


and that they will be skew convertible if 
a+d=0, +d =0, aa +b +b'c+dd'=0, 


the first two of these equations being the conditions in order that the two matrices 
may be respectively periodic of the second order to a factor près. 


45. It may be noticed in passing, that if Z, M are skew convertible matrices of 
the order 2, and if these matrices are also such that Z?=—1, M?=-— 1, then putting 
N=LM=-— ML, we obtain 

P=-1, M?=-1, N?=-1, 
L=MN=-NM, M=NL=-NL, N=LM=-ML, 
Which is a system of relations precisely similar to that in the theory of quaternions. 


62—2 
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46. The integer powers of the matrix 
M=(a, 6b), 
| oad| 
are obtained with great facility from the quadratic equation; thus we have, attending 
first to the positive powers, 
M = (a + d) M — (ad — be), 
M =[(a + d} — (ad — be)] M — (a + d) (ad — be), 
&IC., 


whence also the conditions in order that the matrix may be to a factor près periodic of 
the orders 2, 3, &c. are 


a+d = 0, 
(a + d? — (ad — bc) = 0, 
&e. ; 


and for the negative powers we have 
(ad — be) M+" =— M + (a + d), 
which is equivalent to the ordinary form 
(ad — be) M = d, —b ); 
bbe 


and the other negative powers of M can then be obtained by successive multiplications 
with M~. 


47. The expression for the nth power is however most readily obtained by means 
of a particular algorithm for matrices of the order 2. 


Let h, b, c, J, q be any quantities, and write for shortness R=—h?—4bc; suppose 
also that h’, b’, c’, J’, q’ are any other quantities, such nevertheless that h’: b’: c =h:b:c, 


and write in like manner R’=—h?—4b’c’. Then observing that a x, r 7R are 
respectively equal to -a . sy : i , the matrix 
( h 2bJ ) 
J (cot q- ~=], == 
(° 4 Tp) VR 

2cJ 

Apc 

VR SSA We 


contains only the quantities J, q, which are not the same in both systems; and we 
may therefore represent this matrix by (J, q), and the corresponding matrix with 
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k, bV, d, J’, g by (J, qg). The two matrices are at once seen to be convertible (the 
assumed relations h’ : b : c’=h: 6b: c correspond in fact to the conditions, 


a—d:0:c¢=a-—d:b: ce, 
of convertibility for the ordinary form), and the compound matrix is found to be 
sin (g + q’) as 
karr gh 7) ; 


and in like manner the several convertible matrices (J, q), WJ’, q), (J”, q) &c. give 
the compound matrix 


sin (q+ +g...) 1 T / ” ) 
(aaa a Fan Gee ia i as: 2d 


48. The convertible matrices may be given in the first instance in the ordinary 
form, or we may take these matrices to be 


E A a ES Cb Be. 

| c, d | | c, a | | oy a’ | 
where of course d—a: b:c=d'—-a’: b : ¢=d"—a": P : c’=&. Here writing 
h=d—a, and consequently R=—(d—a)—4be, and assuming also J=4VR and 
cot q “te , and in like manner for the accented letters, the several matrices are 


respectively 


AVR, DGE, q), VR’, q’), &e., 


and the compound matrix is 


sin(g +g +9" ---) a VR)(AVR) (AV R" rarest 
ESSET AR. RYGVE)VGVR)..., q+ +g +). 
49. When the several matrices are each of them equal to 
(a, b ), 
|o a | 


we have of course q=q¢=q'..., R=H=R’..., and we find 
( a, b » sin ng V Ryn é 
|o, a | 7 (amg OYP 4): 


r substituting for the right-hand side, the matrix represented by this notation, and 
Putting for greater simplicity 


SM gWRy= (VE) L or b= VEJ, 


nn 
sin” 7 
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we find ae 
(a, b P=(4L(WReotng—(d—-a)), Lb ) 
| c, d | | Le ; 4 L(~Rcotng+(d—a)) | 
where it will be remembered that 
R=- (d — a} — 4bc and cot g= TE, 
the last of which equations may be replaced by 
—. V—R 
goa Wr Pe guilt 
cos g + sin q Ot = pe 


The formula in fact extends to negative or fractional values of the index n, and when 
n is a fraction, we must, as usual, in order to exhibit the formula in its proper 
generality, write g+2m7 instead of q. In the particular case n= 4, it would be easy 
to show the identity of the value of the square root of the matrix with that before 
obtained by a different process. 


50. The matrix will be to a factor près, periodic of the nth order if only sin ng =0, 


that is, if qa (m must be prime to n, for if it were not, the order of periodicity 
d+a 


would be not n itself, but a submultiple of n); but cos q airs on and the condition 
ad — be 


is therefore 
(d + a)? — 4 (ad — be) cost = 0, 
or as this may also be written, 3 
d? + a? — 2ad cos =m + 4be cos? “zt = 0, 


a result which agrees with those before obtained for the particular values 2 and 3 
of the index of periodicity. 


51. I may remark that the last preceding investigations are intimately connected 


ax+b 


with the investigations of Babbage and others in relation to the function da = E 


I conclude with some remarks upon rectangular matrices. 


52. A matrix such as 


where the number of columns exceeds the number of lines, is said to be a broad 
matrix; a matrix such as 


[Meee 
a. b’ 
a”, 5” 


where the number of lines exceeds the number of columns, is said to be a deep matrix. 
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53. The matrix zero subsists in the present theory, but not the matrix unity. 
Matrices may be added or subtracted when the number of the lines and the number 
of the columns of the one matrix are respectively equal to the number of the lines 
and the number of the columns of the other matrix, and under the like condition 
any number of matrices may be added together. Two matrices may be equal or 
opposite, the one to the other. A matrix may be multiplied by a single quantity, 
giving rise to a matrix of the same form; two matrices so related are similar to 
each other. 


54. The notion of composition applies to rectangular matrices, but it is necessary 
that the number of lines in the second or nearer component matrix should be equal 
o the number of columns in the first or further component matrix; the compound 
matrix will then have as many lines as the first or further component matrix, and 
s many columns as the second or nearer component matrix. 


55. As examples of the composition of rectangular matrices, we have 
( a,b,c § aoea. ) = ( (a, b cha, e’, i), (a, b, ctos f', 7’) (a, b, ote, g, k’), (a, b, ehd, k t) ), 


| d, e,f | ef’; gy; h’ | (d, ef Xa’, e, w) (d, e, f XU, JY 7) (d, e f Xe’, g, k’), (d, e, fF id’, h’ v) 
i, pA k’, j' 


and 
( maA a’, 0, ¢, d) ) =((a, dýa', e), (a, dÝ’, f'), (a, dhe’, 9’), (a, dX’, h’) ). 
b,e | oUF, gt (b, e\a’, e), (b, e Xb, f’), (b, eke’, 9’), (b, eXd’, h’) 
c, f (co, FXa, e), C, FRO, F), (6 FR, g), ( FAY, w) 
56. In the particular case where the lines and columns of the one component 
matrix are respectively equal in number to the columns and lines of the other com- 
ponent matrix, the compound matrix is square, thus we have 
( a, b, ¢ v a’, d )=( (a, b, cha’, b, c’), (a, b, cid’, č, JÀ ) 
[ae f|| v,e | |@e, fX, V, e), (d, e FXX, e, f)| 
c’, £ 


(a, d y a, b, c )=( (a, dYa, d), (a, d'%b, e), (a, d’'Ye, f) 
|d e f| | Œ, e Xa, d), &, e Xb, e} ©, e Xo, f) 
o, f' C, f'Ka dy C, FE, 0, C, fA S) 
The two matrices in the case last considered admit of composition in the two different 
orders of arrangement, but as the resulting square matrices are not of the same order, 


the notion of the convertibility of two matrices does not apply even to the case in 
question. 


and 


57. Since a rectangular matrix cannot be compounded with itself, the notions of 
the inverse or reciprocal matrix and of the powers of the matrix and the whole resulting 
theory of the functions of a matrix, do not apply to rectangular matrices. 
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58. The notion of transposition and the symbol tr. apply to rectangular matrices, 
the effect of a transposition being to convert a broad matrix into a deep one and 
reciprocally. It may be noticed that the symbol tr. may be used for the purpose of 
expressing the law of composition of square or rectangular matrices. Thus treating 
(a, b, c) as a rectangular matrix, or representing it by (a, b, c), we have 


tre( a’ 0, ENS CO), 
| | 
b 


A 


Cc 


and thence 
( a, b, c ) tr. ( a’, WY, ak, b, jji a’ )=(a, b, cha’, b’, c’), 
| | 
b’ 


c 


so that the symbol 
(a, b, cha’, b’, c’) 


would upon principle be replaced by 
( a, b, ¢ ) tr. (a, b’, c ): 
it is however more convenient to retain the symbol 
(a, b, cXa’, b’, c’). 

Hence introducing the symbol tr. only on the left-hand sides, we have 

a, b,c ) tr. (a,b,c )=( (a, b, cha’, vc’), (a, b, cd, e, f) ), 

E Eau, KEANE EGTA A E TA e A 
or to take an example involving square matrices, 


(ab) tr (a,b )=( (a, ba, b), (a, bE, e) ); 
|d, e | | OAA hoy | (d, ea’, b), (d, ed’, e) | 


it thus appears that in the composition of matrices (square or rectangular), when the 
second or nearer component matrix is expressed as a matrix preceded by the symbol 
tr., any line of the compound matrix is obtained by compounding the corresponding 
line of the first or further component matrix successively with the several lines of the 
matrix which preceded by tr. gives the second or nearer component matrix. It is clear 
that the terms ‘symmetrical’ and ‘skew symmetrical’ do not apply to rectangular 
matrices. 
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